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Abstract 

We suggest to use operator product expansion and QCD sum rule techniques to 
estimate the long-distance contribution to the amplitude of the exclusive rare decay 
B -^ K*^. In contrast to a phenomenological description in terms of ip resonances 
converting into a photon, the virtual charm quark loop interacting with soft gluons is 
represented by a new quark-gluon-photon operator. The matrix element of this operator 
is calculated in the same approximation as the matrix element of the leading magnetic 
penguin operator. The overall correction is found to be small, not more than about 5%. 
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1 Introduction 

Rare decays of B mesons such as -B ^ K*'y discovered at CLEO |I|, proceed mainly through 
loops, and are therefore an important tool to study new physics, once the standard model 
contributions are sufficiently understood. The generally accepted calculational procedure is 
to use an effective Hamiltonian 

n^ff = 4%K.v;: y: a(/i)o.(/i) (1) 

with suitable operators Oi and Wilson coefficients Ci, renormalized at the scale /i. The task is 
then to evaluate the Wilson coefficients and the matrix elements of the operators between the 
initial and final hadronic states, in the case at hand, between B and K*. The coefficients C, can 
be calculated perturbatively. A comprehensive review of the current status of the Hamiltonian 
(|I|) can be found in 0. In order to obtain the exclusive matrix elements {K* \ Oi \ B), some 
hadronic model or nonperturbative method must be invoked. 

In leading logarithmic approximation, only the matrix element of the dominant O7 operator 
is needed. It has been calculated using various methods [^ §, ^, ||, |^. However, there remains 
a rather strong /i dependence. It can be reduced by next-to-leading order calculations of 
the Ci and a corresponding higher order in Og determination of the matrix elements. For 
example, in the case of the inclusive reaction b ^ s'j one has to calculate the two-loop- 
diagrams involving the operator O2. Here, the next-to- leading analysis has essentially been 
completed [^, |^, ^ |TT|. Generally, it is assumed that a similar approach is adequate for 



exclusive transitions. The corresponding diagrams are shown in Fig. 1. 

When considering the penguin-like loops in these diagrams, a perturbative calculation 
seems adequate for gluon virtualities j/c^l > ml and even when ml < \k'^\ < ml. Furthermore, 
if the photon is emitted from the 6 or s quark or from the spectator quark (Figs. lb,d), 
the integration over the c-quark loop produces a factor k"^ which cancels the denominator 
of the gluon propagator, effectively suppressing low-virtuality gluons. However, when the c- 
quark loop converts into a photon (Fig. la,c), the gluon may have a characteristic virtuality 
|fc^| -C Am^l. This case is beyond perturbation theory. The role of these long-distance penguins 
in B ^ K*'-f illustrated graphically in Fig. 2 will be the main concern of this paper. 

Their effect has been considered by many authors on a phenomenological level. The pos- 
sibility extensively discussed in the literature [|l^] is to model the transition B — ^ M7, M 
being a light meson, by a sum over the nonleptonic decays B — »• V*M followed by the transi- 
tion V* —>■ J where V* is an appropriate virtual vector state such as the J/ip. Since at zero 
momentum squared p^ of the photon, the characteristic distance to the physical threshold 
p^ = M| is rather large, not only the ground state J/ip but also excited ip states must be 
included. Unfortunately, the transverse parts of these amplitudes are poorly known and the 
extrapolation from p^ = Ml to p^ = is not straightforward. Therefore, for a real photon, the 
vector-resonance description of the virtual cc loop seems uncertain, calling for a more reliable 
or at least different estimate. 

As a new approach to the problem, we suggest to use operator product expansion (OPE) 
and QCD sum rules [0]. Several important circumstances make these methods exception- 
ally suitable. Firstly, the fact that the cc loop is far off-shell, which is a disadvantage for ip 
insertion, turns out to be a major advantage here. It means that we must keep only a few 
terms in the OPE. Secondly, the QCD sum rule method provides a physically very intuitive 
parametrization of the interaction of the cc loop with soft gluons in terms of vacuum conden- 
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sates. Thirdly, and most importantly, we will be able to compare the long-distance effect in 
B — *■ K*^ with the contribution of the dominant magnetic penguin operator obtained in one 
and the same calculational framework. Of course, as with any low energy effective method, 
there are some unavoidable theoretical uncertainties. One of them is the question of double 
counting in the perturbative and nonperturbative contributions. As usual, we assume that 
the former are small in the region where the latter dominate, that is at \k'^\ <^ 4m^. 

There are other nonperturbative effects, such as those illustrated in Fig. 3, which involve 
radiative matrix elements of the four-fermion operators Oi2- These contribute already at tree 



level and have been estimated ||IJ, |T^ using QCD sum rules on the light-cone. They are 
unimportant for the b ^ s transitions due to CKM suppression (Fig. 3a), and also relatively 
small for b ^ d transitions (Fig. 3b). 

The paper is organized as follows. In Sect. 2, we consider the contribution of the operator 
O2 to the exclusive B -^ K*^ amplitude and single out the term responsible for soft gluon 
emission from the c-quark loop. With the help of OPE, this term is expressed by a composite 
operator involving quark, gluon and photon fields. In Sect. 3, we derive a QCD sum rule for the 
B -^ K*^ matrix element, including the new operator. The most important nonperturbative 
contributions of gluon, quark-gluon, and quark condensates up to dimension 6 to the relevant 
correlation function are calculated in Sect. 4. Our numerical results are presented in Sect. 5 
and our conclusions in Sect. 6. 

2 Long-distance contribution of the operator O2 

In the framework of the effective Hamiltonian (|I|) the exclusive B -^ K*'y decay amplitude is 
given by 

AiB -. K*7) = {K*{q) 7(p) I K// I Bip + q)) = 



AG 



F 



- V,tV:,{K*{q) I [C,Oj + ze>^ ^ C J d'x e^^^T{j7(x)O.(0)}] | Bip + q)) , (2) 

where the four-momenta of the B and K* meson, and the photon are denoted hj p + q, q, and 
p, respectively, while e^ is the photon polarization four-vector. The electromagnetic current 
is given by 

Jr= E e,g7,g , (3) 

q=u,d,s,c,b 

where e^ are the quark electric charges. 

In (0), we have isolated the operator O7, as it is the only operator which contains the 
photon field at the tree level: 

O7 = ^.-SLcr^^bRF^" , (4) 



with qjiL = |(1 ± 75)g, F^^, being the photon field tensor, and e = yiyra being the electro- 
magnetic coupling. Note that throughout the paper we put m^ = 0. The other operators 
generate radiative decays when combined with the electromagnetic interaction of quarks with 
the photon. In this paper, we restrict ourselves to the four-quark operators 

Oi = (cl7a.cl)(sl7^&l) , (5) 

O2 = {sLl^.CL){cLl^bL) , (6) 



relying on the fact that the other operators have much smaller short- distance coefficients. To 
our knowledge, an analysis of long-distance effects in B ^ K*'j generated by O^^q and Og 
has never been done and is certainly desirable. Perturbative contributions of Og have been 
estimated in the relativistic quark model |T^ and found to be small. 

For O2, the relevant diagrams contributing to the amplitude (§) are shown in Fig. 1. The 
analogous graphs for the operator Oi vanish for a real photon. After Fierz transformation, 
the operator O2 reads: 

02 = lo, + Id, , (7) 
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di = {cT^—c){sT^—h) , 



where F^ = "^^{1 — 75). It is only the operator Oi which plays a role. The effective colour 
transfer in the weak vertex is compensated by a gluon exchange between the c-quark loop and 
other quarks in the process. The photon emission from the c-quark loop is therefore a purely 
nonfactorizable effect. 

Keeping in (|^) only the operators Oj and Oi and retaining in j^J^ only the c-quark part, 
the decay amplitude simplifies to 



A{B -. K*^) = ^V,,V:,{K*{q) 



C7O 



■jKj-j 



2 

A 



+ i^e^' I d^x e'P"T{c(a;)7^c(x)di(0)} | B{p + q)). (9) 



In this expression, one encounters a product of C7^c and cFp— c currents generating a c-quark 
loop. When the photon is on-shell, the c-quark pair is far off-shell. In this situation, we expect 
that the OPE for the c-quark loop in terms of local operators is adequate. 

The ffist term of this expansion is simply a gluon field operator. The corresponding Wilson 
coefficient is represented by the diagrams of Fig. 4. It can easily be calculated using the fixed- 
point {x = 0) gauge which in the approximation needed reads 

Al^ix) = l/2x^G%{0) . (10) 

The OPE then results in 

e^^^T{c(x)7;.c(x), c(0)F,yc(0)} = C^p.^gG'^'''' + ... (11) 

with the short-distance coefficient {p^ = 0), 
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CalBt^p = - 2 2 (P«Pp^/3m + PaPt^gfip) , (12) 

t:0 /I 1 1 L ^ 

and G°-°'f^ = l/2e°^^^G^;^. Note that in the diagrams of Fig. 4 only the axial part of FpA" 
contributes, whereas the contribution of the vector part vanishes due to Furry's theorem. The 
higher dimension operators in the expansion (|TT]) denoted by the ellipses contain derivatives of 
the gluon field and/or additional gluon field operators. They are neglected here. Substituting 
the expansion (^) into (^, and putting Cc = (2/3)e yields 

A{B ^ K*^) = 4 (^] V,bV;AK*{q) \ [G,Oj + G^Op] \ B{p + q)) (13) 



where 



O, 



2887r2m; 



:sT,-gGl0DPF^^b 



(14) 



is a new effective operator describing the soft gluon interaction sketched in Fig. 2. 

Before turning to the estimate of the matrix element of Op, a comment on the accuracy 
of approximation (|TT]) is in order. A similar approximation was used in |T^ in estimating 
the nonfactorizable contributions to i? — ;> Dn (see also |l^). In this analysis, the criteria for 
the applicability of OPE are discussed in detail. Following the same line of arguments, we 
notice that the next-to-leading operator in ( pTD will contain a derivative of the gluon field DaG 
multiplied by the photon momentum pa and an extra factor of order of l/4m^. The overall 
suppression factor will thus be of order of yU/iPo/4m^ where fit is a characteristic hadronic 
scale, and po ^ fnii/2 is the photon energy. Assuming /i^ ~ 300 MeV, mj, ~ 5 GeV and 
rric — 1.3 GeV one expects a 10% effect. From this very rough estimate the truncation of ( [TT| ) 
after the first term appears justifiable. On the other hand, one clearly cannot apply the same 
approximation to light-quark loops. 

Next, we write the matrix element of the magnetic penguin operator in the form 



{K*{q)\Oj\B{p + q)) 



enib 
327r2 



[^^aPc 



.g^-2(7.„g^)e^,F-^. 



(15) 



where ck* is the polarization vector of the K* and F^,^ = i{eyp^ — e^py). In the above, the 
relation 



C/xi^TS — —-^tivapO' 



aP 



(16) 



has been used, which leaves only one independent invariant amplitude F. Similarly to (|TBp 
we parametrize the matrix element (|13D of the operator Op in the form 



{K*{q)\OF\B{p + q)) 



2887r2m2 



{iLeaparq' - 2Lg^aqi3) e^.F"'' 



(17) 



Note that the invariant amplitudes L and L have dimension GeV^ and that, in general, L ^ L. 
Finally, substituting (|T^) and (|T3p in (|TB|) one gets 



A{B -^ K*-f) 



'eGFVcbV:, 



-K" 



iea/sarq'^ I -GjrribF + -^r^L 



+ '^gac.qp 



-G-rrUbF + — 4l 



9ml 






-K 



While the amplitude proportional to F has been investigated previously, the long-distance 
corrections proportional to L and L are new. 



3 QCD sum rule for the B -^ K*^ amplitude 

Estimates for the matrix element (|I5|) already exist in the framework of local 0, ^ and light- 
cone QCD sum rules . The additional matrix element (|17D can in principle be calculated in 
the same framework with the same accuracy. The light-cone approach to this matrix element 
involves quark-antiquark-gluon wave functions of the K* meson. Since these wave functions 



are not yet sufficiently known, we shall estimate ([T7|) using the local OPE. Following [^ |^, 
we consider the three-point correlation function 

T:^(p,g)F°^ = -^Jd^xd%exp[-z{p + q)x + tqy]{0\T{u{y)Ys{y), 

[C7m,.-(0)a„^(l + 75)&(0) - ^Pp3(0)r''^^G^^(0)6(0)], 6(a;)z75«(x)}|0)F"^ (19) 

interpolating the complete matrix element in (0), with the B and K* mesons replaced by 
appropriate generating currents. Moreover, for the operators O7 and Op we have substituted 
the explicit expressions (H) and ([T^) . Isolating in (|^) the tensor structures (^) and ([T7|) we 
write 

T^gF'^^ = ^ {-Tie^p^rQ' F^P + 2fqpF,p + ...) , (20) 

where the invariant amplitudes T and T are functions of the two independent variables {p + qY 
and q^. 

Inserting in the r.h.s. of ([19|) the complete sets of hadronic states carrying the B and K* 
quantum numbers one obtains the double dispersion relation 

T{{ + f 2n m^mK^fsfK* (c m F ^^ l\ 

mfe(m| - (p + g)2)(m|^, - g2) \^ 9^2 J 

+ / dsids2- ^ (21) 

J [si - [p + qY)[s2 - q^) 

where we have used (^5]) and ([T7|) together with 



2 
{B I 6^75M I 0) = ^/b , (22) 



(0 I u-fyS I /C*) = mK*fK*eK*u ■ (23) 



The integral over the double spectral density ph in ( PT| ) represents the contributions of excited 
S and K* states as well as of continuum contributions. Possible subtraction terms are omitted 
for brevity. The analogous dispersion relation for T is obtained from (^) by the replacements 
L ^ L and ph —> ph ■ 

Furthermore, at p^ = and at large spacelike values of {p + qY and g^, the amplitudes 
T and T can be evaluated in QCD in terms of vacuum expectation values of certain local 
operators VL^. 



TQCD{{p + qY,q')=T. 



Crm.Uip + g)^ q') - ^V,{{p + q)\ q^) 



W . (24) 



A similar representation holds for Tqcd with Ud replaced by f/^. The index d refers to the 
dimension of the operators. In perturbation theory, only the d = unit operator contributes 
to (P^. The nonperturbative effects are accumulated in terms with d > 0. In this analysis, 
we take into account (^3) = (gg), (^4) = {^C'^^G"^") , (^5) = mHqq), and (fie) = {qgY- 
These are the standard parametrizations for gluon and quark-gluon condensate densities. For 



the four-quark condensate the vacuum insertion approximation [13| is used 



Following the usual procedure, one equates (pTD and (0), and replaces the double spectral 
density p^ in (^) by the double imaginary part of the amplitude Tqqd- Finally, the Borel 
transformation 



BM^fiQ^)^^^ lim 



{Q 



2\n+l 



d 



n\ 



dQ^ 



fm ^ f{M' 



(25) 



with respect to the variables {p + q)^ and g^ is applied, in order to suppress exponentially the 
contribution of higher states and to get rid of subtraction terms in the dispersion relations. 
The second invariant amplitude T is treated in the same way. 
This yields the sum rule 



CrTTibF 



a, 

9m^ 



:L 



rrib 



m%mK*fBfK* 



^ml/Ml+m\JMl 



X' 



vr^ 



K^ 



•} 



dsids2 



+ E 



C^mJmsJms^To^Si, S2) 



C2 
9m? 



ImsJms,^Uo{si, S2) 



C-rmMMl Ml) - ^Ud{Ml M, 



W . 



-si/Aq-S2/Ml 



(26) 



where the threshold parameters Sq and s^* indicate the subtraction of higher state contribu- 
tions. A second sum rule is obtained from (^) by replacing L ^ L and Ud -^ Ud- 

In the approximation considered, the above sum rule can obviously be divided into two 
separate sum rules, one for F (terms involving T^ ) and one for L (terms involving Ud)- The 
sum rule for F have already been analyzed ||^ including terms up to d = 6, in lowest order in 
as- In the following, we calculate the sum rule for the long-distance corrections L and L in 
the same approximation. 



4 Calculation of the long-distance amplitudes 

We need to calculate the coefficients Ud and f/^ up to (i = 6. The coefficients f/o, f^o, U?^ and 
f/3 are of order Ug as is obvious from Fig. 5. Since the 0{as) corrections to the sum rule for F 
(exemplified in Fig. 5a, d) are not taken into account in we put f/o = [/q = ^/s = f/3 = in 
(^61) . A systematic inclusion of these effects would require two-loop calculations and a proper 
treatment of operator mixing. Thus, to the requested accuracy, we have 



rrib 



ni%mK*fBfK 



_^ml/Ml+m\JMl 



[u,{MlMl)(^Gl,G^n 



+ U,{MlMl)ml{qq) + U,{MlMl){qqf] , 



(27) 



and a similar expression for the amplitude L with f/4,5,6 replaced by f/4,5,6- 

In order to calculate the Wilson coefficients t/4 and f/4, one has to contract all quark 
fields in the second term of the correlation function (|19|). One of the three resulting quark 
propagators is substituted by the first order expression in the external gluon field AJ^, while 
for the other two propagators the free approximation is used. This yields the diagrams of 
Fig. 6a, b, c. The calculation of these diagrams is performed in the Fock-Schwinger gauge 
(pni). Furthermore, when expressing the result in form of a double dispersion integral in the 



variables {p + q)^ and g^ one may omit all terms which depend only on one of these variables 
since such terms vanish after Borel transformation. The double dispersion contribution of the 
diagram Fig. 6c vanishes altogether. The result is 



U, 



jj{Fig.6a) jj{Fig.6b) 



(28) 



with 



jj(Fig.6a) 



'48 



dsi 



(p+qy 



ds-} 



U. 



(Fig.Gb) 



48 



ds^ 



Si - (p + qy 



S2 -q^ 

dS2 



nit 



Sl -S2 



S2 -g' 



+ 



and 



U, 



{S2 -q^){si - S2) 



jj{Fig.6a) jj{Fig.6b) 



1 



+ 



Sl 
1 



S2 



Sl - (p + g)2 S2 - g 



, (29) 

(30) 

(31) 

5 quark-gluon condensate term are obtained from 
Fig. 6d which involves the vacuum average of the uuG operator. The analogous diagram, 
involving ssG depends only on g^ and can therefore be dropped. One finds 



with 



The coefficients f/5 and U^ of the d 



{Fig. 6a) _ jj{Fig.6a) jj{Fig.6b) ^ ^{Fig.Gb) 



U, = U, 



ml 



12(m2-(p + g)2)(-g2) • 



(32) 



There are two d = 6 contributions to the correlation function ([191) . One of them originates 
from expanding u{y) and u{x) at x = y = 0. In vacuum insertion approximation, the vacuum 
average of the resulting operator uVpuGI^^ is expressed in terms of (gg)^ according to the 
1l9| . The corresponding coefficients derived from Fig. 6d read 



formula given in 



U, 



(Fig.Gd) 
6 



f/. 



(Fig.ed) 



AnasTTih 



ml/q^ 



27 (-g2)(m2-(p + g)2) 



(33) 



The second d = Q contribution comes from the four-quark condensate. The relevant diagrams 
are shown in Fig. 6e and 6f. The coefficients are 



f/. 



(FigSeJ) 







^(Fig.GeJ) 



Airasmi, 1 — ml/q^ 



9 



q^){ml - {p + qy 



In total, we have 



U. = U, 



{Fig.Qd) 



(34) 

(35) 
(36) 

After double Borel transformation (^) of the above expressions for Ud and Ud in {p + qY 
and g^, and substitution in (|27| ) and in the analogous expression for L, we obtain the final 
sum rules 



U. = U, 



{Fig.Gd) 



+ u, 



(Fig.GeJ) 



rUb 



mK*m%fBfK^ 



exp 



m 



B 



Mf 



+ 



ni 



K* 



Ml 




(irexp(— r) 



rmr 



T 



Mf 



Ml 



Ml 



1- 



mr 



rMf 



1 + 



Ml 

Vm! 



ml{qq)ml A-nas{qq)'^mb 



12 





mK*m%fBfK 
drexp 

ml{qq)ml 16nas{qqymb ( m 



12 



27 



1 + 



Ml 



exp 



ml 

'Ml 



(37) 



(38) 



5 Numerical estimates 

For the numerical estimate of the amphtudes L and L we take the values of the physical 
parameters also used in 0, that is m,, = 4.8 GeV, {^GG) = 0.012 GeV^ {qq) = -(240MeV)3 



at the scale /i 



1 GeV, ml 



0.8 GeV^, A 



(4) 
QCD 



230 MeV, Jk* = 210 MeV, and fe = 140 



MeV. The latter value results from the corresponding two-point sum rule without 0{as) 
corrections and should not be confused with the physical value of fs which is larger. We also 
use the same interval of the Borel mass Mi, 7 < Mf < 10 GeV^, keeping the ratio Mf/M^ = 3 
fixed. The characteristic scale at which the amplitudes are estimated by the sum rules is of 
order of the average Borel parameter that is, of order of the virtuality of the quarks in (plOl) . 



In 0, the scale /is = (Mi Ms 



,1/4 



2 GeV is chosen for evaluation of F. Neglecting an 



inessential evolution of the quark-gluon condensate and of the product dsiqq)^ from fJ, = I 
GeV to fi = fiB, the sum rules (^) and ( PS]) yield: 



(39) 



L = 0.55 ±0.1 GeV^ 



L = 0.70 ±0.1 GeV^ 



The theoretical uncertainty quoted here reflects only the variation of L and L with the 
Borel masses Mi and M2 in the accepted interval. In contrast to 0, we have neglected 
the finite strange quark mass effect which however makes no difference numerically. The 
dominant contributions in ( pTj ) and ( p^ ) come from the quark-gluon condensate (75 %) and 
gluon condensate (15-20 %), while the d = 6 terms do not exceed 10%. The deviation of L 
from L which is numerically not very significant is generated partly by the gluon condensate 
and partly by the four-quark condensate. 

For comparison, we quote the corresponding prediction for the penguin amplitude (at 
fi = fis) given in 0: 

F = 0.9 ±0.1. (40) 

With (|39D and (0) we find for the relative magnitude of the long-distance corrections in (plsl) : 



C, 



L 



9C7/ mi,mlF 



0.02^0.04 



c 

^ ' ^ 0.03^0.05, (41) 



,— 9C7/ nihrnlF 

Apart from rrih = 4.8 GeV and rric = 1.3 GeV we have used here C2 = 1.14 and Cj = —0.30, 
that is the leading-logarithmic approximation of the short-distance coefficients at /i ~ mb 



Since both estimates ( pQ]) and (^OD are obtained by the same method and in the same 



approximation, one can hope that many of the uncertainties cancel in the above ratios. 

6 Conclusion 

In this paper, we have estimated the long-distance effect in B —^ K*'y due to interactions of 
the virtual c-quark loop with soft gluons. Using OPE we have derived a new quark-gluon- 
photon operator Op given in (|l^ which is responsible for the effect. Furthermore, we have 
obtained a QCD sum rule for the matrix element {K*'y \ Op \ B) taking into account the most 
important nonperturbative contributions. Quantitatively, the correction is found to increase 
the magnetic penguin amplitude by a few percent. This result indicates that long-distance 
effects do not play an important role in the exclusive decay B — > K*'-). 

The estimate presented here is not free from theoretical uncertainties. They are partly 
due to the lack of knowledge of the perturbative 0{as) corrections to the Wilson coefficients. 
Related to that is the ambiguity in the choice of the scale /i. In addition, there are uncertain- 
ties inherent to the short-distance sum rules. It would therefore be desirable to repeat the 
calculation using light-cone sum rules for a cross-check. However, the latter approach requires 
a better understanding of the higher twist light-cone wave functions of the K* meson. 

When this work was completed, the paper |20| appeared in which the nonperturbative 



correction to the inclusive B -^ Xgj rate due to soft gluon emission by the c-quark loop is 
estimated. The responsible effective operator obtained in ^^ coincides with our operator 



|T^. After revision of pO|, the effect predicted for the inclusive rate is similar in size to the 



correction found here for the exclusive channel B -^ K*'~f, however, opposite in sign. 
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(d) 



Figure 1: Diagrams contributing to the B -^ K*'~f amplitude generated by the operator O2 
(black circle). 




Figure 2: Interaction of the virtual c-quark loop with soft gluons. 




(a) (b) 

Figure 3: Photon emission by quarks at long-distances in (a) B -^ K*'j and (b) B ^ fry. 
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Figure 4: Diagrams contributing to the operator product expansion of the c-quark loop. 






(a) 
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(e) 



(f) 



Figure 5: Diagrams determining the 0(0;^) corrections to the Wilson coefficients of the three- 
point correlation function (|l^). The sb vertex involves O7, while the sbG vertex involves 
Of. 
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(d) 
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Figure 6: Diagrams yielding the coefficients (a-c) f/4, f/4 , (d) f/5, U^, and (d-f) Uq, Uq in the 
sum rules for L and L. 



12 



